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^ ' Analytic solutions are presented which describe the construction of a traversable 

o' 

^' ■ wormhole from a Schwarzschild black hole, and the enlargement of such a wormhole, 

0^ I in Einstein gravity. The matter model is pure radiation which may have negative en- 

ergy density (phantom or ghost radiation) and the idealization of impulsive radiation 
(infinitesimally thin null shells) is employed. 
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^ : I. INTRODUCTION 



While black holes are now almost universally accepted as astrophysical realities, 
j^ ■ traversable wormholes are still a theoretical idea [l|, \^. Yet they are both predictions of 

General Relativity in a sense, though black holes require positive-energy matter (or vac- 
uum) whereas wormholes require negative-energy matter. While normal positive-energy 
matter was long thought to dominate the universe, it is now known that this is not so. The 
recently d.scove.ed acceleration of the un.ve.se SQ ind.cates that its evo.ut.on is dominated 
by unknown dark energy which violates at least the strong energy condition (w > — 1/3 to 
cosmologists, where w is the ratio of pressure to density in relativistic units, for a homoge- 
neous isotropic cosmos), and perhaps also the weak energy condition {w > —1), where it is 

n 

known as phantom energy ^5|. Such phantom energy is precisely what is needed to support 
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traversable wormholes (a lD, |8|, |9| . 

While black holes and traversable wormholes have been regarded by most experts as quite 
different, one of the authors has argued that they form a continuum and are theoretically 
interconvertible ^. Specifically, both are locally characterized by trapping horizons a, 
y, llO, lll|, which are the Killing horizons of a stationary black hole and the throat of a 
stationary wormhole. The difference is the causal nature, being spatial or null for a black 
hole and temporal for a wormhole. This in turn depends on whether the energy density 
is positive, zero or negative. If the energy density can be controlled, it should be possible 
to dynamically create a traversable wormhole from a black hole and vice versa. This was 
first concretely demonstrated in a two-dimensional model |12] , but it is more difficult in full 
General Relativity. Numerical simulations have been used to study a wormhole supported 
by a ghost (or phantom) scalar field, showing that it does indeed collapse to a black hole 
if perturbed by positive-energy matter (l3|. (We use phantom to mean that the energy 
density has the opposite sign to normal, which is equivalent in our cases to the convention 
for ghost fields in quantum field theory, that the kinetic energy has the opposite sign). 
As for analytic results, one simple case is a static wormhole supported by pure ghost (or 
phantom) radiation 14] ; it is easy to see that if the radiation is switched off, it immediately 
collapses to a Schwarzschild black hole. The converse, creating a traversable wormhole from 
a Schwarzschild black hole, is more complex and is the first main result of this article. 

As above, we use pure phantom radiation as the exotic matter model. We also employ 
the idealization of impulsive radiation, where the radiation forms an infinitesimally thin 
null shell, thereby delivering finite energy-momentum in an instant jl^. Space-time regions 
can be matched across such shells using the Barrabes- Israel formalism [16[. This allows an 
ingenious analytic construction of the desired type of solution, by matching Schwarzschild, 
static- wormhole and Vaidya regions, the latter consisting of pure radiation propagating in a 
fixed direction |l2|. Our other main result is the similar construction of analytic solutions 
describing the enlargement or reduction of such a wormhole. Then if Wheeler's space-time 
foam picture [1^ is correct and Planck-sized virtual black holes are continually forming, we 
have exact solutions in standard Einstein gravity describing how they may be converted into 
traversable wormholes and enlarged to usable size. Our results have been summarized in a 
short, a*. H. 

Our paper is organized as follows. In Sec. m we briefly review the static wormhole 



solutions with pure phantom radiation, and the Schwarzschild and Vaidya solutions. In 
Sec. mil we show how to join these basic solutions at null boundaries. In Sees. |TVl and El 
we present analytic solutions which describe respectively the construction and enlargement 
of a wormhole. In Sec. IVII we consider the jump in energy due to impulsive radiation, as a 
check that the matchings are physically reasonable, and as a simple way to understand the 
changes in area of the wormhole throat or black-hole horizon. The final section is devoted 
to summary. 

II. BASIC SOLUTIONS 

In this section we review the traversable wormhole solution |l^, the Schwarzschild solu- 
tion and the Vaidya solutions in the various coordinates needed. We will consider spherically 
symmetric space-times only. It is convenient to use the area radius r = ^jAj^^ where A 
is the area of the spheres of symmetry. Although we often use r as a coordinate, it is a 
geometrical invariant of the metric and is assumed throughout to be continuous. A useful 



quantity is the local gravitational mass-energy 
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m 



E='^{l-9''\,r,^). (1) 

Note that E = r/2 on a trapping horizon, where g^^r^^r^y = 0, including both black-hole 
horizons |10] and wormhole mouths 8|. The energy-momentum tensor of pure radiation (or 
null dust) is Tat = pUaUb, where Ua is null and p is the energy density. Normally p > 0, but 
p < defines pure phantom radiation. 



A. Static wormhole solution with pure phantom radiation 



Q 



The static wormhole solutions jlj| supported by opposing streams of pure phantom ra- 
diation can be written as 
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ds = -TT21W ^ ^0V^^' + ' ^^ ^'^ 

where t is the static time coordinate and dQ"^ refers to the unit sphere. Here / is a function 
of r, 

r = a(e-''+2Z0), (3) 



and is an error function, 

0(/) = f e-^'di + &, (4) 

Jo 

where a > 0, b and A > are constants. 

The local energy evaluates as E = e where 

e=^(e-''+2/0-2e'V')- (5) 

Using the mass-energy e, the metric Q is rewritten as 

In the case 6 = 0, the solutions (j21) or (jH)) describe symmetric wormholes. The spacetime 
is not asymptotically flat, but otherwise constitutes a Morris-Thorne wormhole. The 6 7^ 
cases include asymmetric wormholes which are analogous to the asymmetric Ellis wormhole 
for a phantom Klein-Gordon field J2j|. For the space-time solution to be a wormhole, the 
inequality 



\b\ <bcr = ^ (7) 

is needed. For any other value of b, a singularity is present. Hereafter, we consider 6 = 0, 
describing a symmetric wormhole with minimal surfaces at the wormhole throat / = 0, with 
area radius r = a. 

The solutions (jHl) niay be written in dual-null form 

2A 

l + 2/e''0 ^ ^ 

where the null coordinates x^ are defined by 

dx"^ = dt±^ (e-^^ + 2l(p) dl. (9) 

Then the radial null geodesies are given by constant x"^. We need the metric (jHI) in both 
ingoing and outgoing radiation coordinates: 

.1^ 






l^fXe 



[\ + 2W\ 



-dx =F 2dr 



+ r'dVi'. (10) 



In these coordinates the radial null geodesies are the lines of constant x^ (choosing one) and 
the curves given by 



dr , a/Ac 



±7T^H^7^- (11) 



dx^ (1 + 2/e'V) 





FIG. 1: (Color online). Penrose diagrams of (i) a Schwarzschild black hole and (ii) a Hayward 
traversable wormhole [lj|. The bold magenta and green lines represent the trapping horizons, 
9+r = and 5_r = 0, respectively, which constitute the event horizons of the black hole and the 
throat of the wormhole. Yellow (light) and gray (dark) quadrants represent past trapped and future 
trapped regions, respectively. Wavy cyan lines represent the constant-profile radiation supporting 
the wormhole structure. 

The Penrose diagram is shown by FigQ (ii). The energy- momentum tensor supporting the 
wormhole is found to be 



T. 



A 



ab 



87ir' 



istK + s-s^). 



(12) 



This is the energy tensor of two opposing streams of pure phantom radiation, with A 
—Airr'^Tu being the resulting negative linear energy density. On the other hand, the solutions 
with pure radiation of the usual positive energy density were founded by Gergely 



B. Sch^varzschild solution 



The Schwarzschild metric is given by 

2M 



ds' 



1 



dr+\i- 



2M 



dr^ + r^dn^ 



(13) 



where the constant M is the Schwarzschild mass, which coincides with the local energy, 
E = M. Rewriting in Eddington-Finkelstein coordinates 



V = t-C{r + 2M ln(l - r/2M)) 



one finds 



ds' = -dV 



2M , 
1 ]dV + 2Cdr 



+ r'dQ' 



(14) 



(15) 
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where C is a sign factor: C is 1 for outgoing radiation or —1 for ingoing radiation, where 
this means that the area respectively increases or decreases afong the future-null generators. 
The Penrose diagram is shown by Fig^ (i)- 

C. Vaidya solutions 

The metric of the Vaidya solutions is given by 



ds^ = -dV 



(l_MlQ),V^ + 2C* 



+ r^dQ\ (16) 



where ^ is a sign factor, where ^ = 1 for outgoing radiation, and ^ = — 1 for ingoing 
radiation. The mass function m coincides with the local energy, E = m. The corresponding 
energy-momentum tensor is given by 

III. MATCHING VAIDYA REGIONS TO STATIC- WORMHOLE AND 

SCHWARZSCHILD REGIONS 

In this section, we derive the matching formulas between Schwarzschild, Vaidya and 
static-wormhole regions along null hypersurfaces, following the Barrabes-Israel formalism 
|lq . This is a preliminary to constructing the wormhole-construction and wormhole- 
enlargement models. 



A. Matching Vaidya and Schwarzschild regions 

Firstly we consider the matching between Schwarzschild and Vaidya regions. We start 
by writing Schwarzschild and Vaidya solutions in the form 

ds^ = -e^dVife'f'dV + 2Cdr) + r^dQ^, (18) 

where the metric functions are 

/5 = 1 - — , ^s = (19) 



for Schwarzschild, and 



;, = 1 _ M10_ „^ _ „ (20) 
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for Vaidya. 

Now we consider the boundary surface V = Vq (constant). The normal to the hypersurface 
^ = V — Vq = is n^ = (a'^dfj,^ = Ca~^6^ , where a is a positive function. (Barrabes 
& Israel took a < 0). For a null hypersurface, the normal is also tangent, so to obtain 
extrinsic curvature one needs a different vector. From n^ Barrabes and Israel introduced 
a so-called transverse null vector A*"^ by requiring N^N'^ = 0, and N^n^^ = —1. Without 
loss of generality, we assume that A''^ takes the form A^^ = NyS^ + NrSJ^, and choose the 
arbitrary function a as a = e~^. Then N^ is given by 

N, = C^< + S;. (21) 

Choosing the coordinates r,6, and if as the three intrinsic coordinates ^" = {r,6,ip), (a = 
1, 2, 3) on the hypersurface \^ = Vq, we find 



e 



(1) ~ "r5 6(2) ~ "e' ^(3) ~ "</" y^^) 



where e'f ■, = dx^ IdE"'. Then, it can be shown that the transverse extrinsic curvature, defined 

7^a. = -N,e\,) {V^e^^ , (23) 

takes the form 

TZa, = diag. |c^, -C y, -C y sin^ e\ . (24) 

The jump in transverse extrinsic curvature is denoted by 

lab = 2[nab]. (25) 



Once lab is given, using the formula 







^ab ^ _gab ^ ^L^g-c^^l^ + gM^a^c _ ^ab^c^d _ gf^b)^^^^ (26) 

iOTT 

we can calculate the surface energy-momentum tensor r'*'' on the null hypersurface V = Vq. 
In components, 

T^>' = atl^ + Pgf, (27) 

where 

gf = r-2 {6X + sin-2 e6X) , 

r = 5^, 1% = 0. (28) 
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Here a represents the surface energy density of the null shell and P the pressure in the Q 
and v? directions. Then 



. = cJa = ,.c''-'"<'°' 



57rr 



47rr^ 



P = - 



c 



^TT 



Or 



0. 



(29) 



We take the sign factor rji to be 1 if the radiation is to the future of the Schwarzschild 
region, and —1 if the radiation is to the past of the Schwarzschild region (Fig|21). The 
energy-momentum tensor of the impulsive radiation is given generally by [l^ 



Tf"" = aT^"'6{<^), 



where 6 denotes the Dirac delta distribution. In our case it reduces to 

rjirr 



r"S(V-V,)=, J":rT"^^ S(V-V,). 



Airr^ 



(30) 



(31) 




.Schwarzschild 



Vaidya 



Vaidya 



Vaidya 




Vaidya 



FIG. 2: Matching domains: iji = 1 (left) and r/i = — 1 (right). 



B. Matching Vaidya and static-wormhole regions 



Secondly we consider matching Vaidya and static-wormhole regions. We rewrite the static 
wormhole solution (fTUI) as 



2jn2 



ds' = -e^duUwe^'^du - 2Cdr) + r'dil 
where u is x^ for ^ = ±1. The metric functions fw and ipw are defined as 



/i 



2e 



w 






(32) 



(33) 
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Now we consider the boundary surface u = Uq (constant). Then, from the previous subsec- 
tion, the transverse extrinsic curvature takes the form 

-?:>+ ^- / r-^^W /-^fw r-^fw ■ 2/3! fOA\ 

Kb = diag. < -C-Q^, C^-, C^- sm ^ ^ . (34) 



On the other hand, we write the Vaidya solution [1T\ in the same form as (^0]) 

ds^ = -e^^dVifve^'^dV + 2Cdr) + r^dQ^, (35) 

where the metric fy and tpy is defined by ()20|) . The hypersurface u = uq in the (V, r) 
coordinates can be written as $^ = \^ — Vo(r) = 0, where Vo(r) is a solution of the equation 

^ = -Cfe-^^ iu = Uo). (36) 

dr fv 

Then the normal to the surface is given by 

n- = Cp-%<^- = CP-' (^6l + |;e-'^"5;) , (37) 

where /3 is a negative and otherwise arbitrary function. From n~ we can also introduce the 
transverse null vector N~ by requiring N^N^^ = 0, and N~n~^ = —1. It can be shown 
that it takes the form 

The basis vectors are 

ei,^ = -CJ^e-^^Si^ + S^, e-^^=8l e"^ = 5^, (39) 

where eV^ = dx't/d^"'. To be sure that the two transverse vectors A*"^ defined in the two 
faces of the hypersurface u = Uq represent the same vector, we need to impose the condition 



^A *2{a) \u=uo — ^X ^(a) \u=uoj (40) 

which requires that the function (3 has to he (3 = — exp{— ?/'_}. Once A^;^ and eVj^ are given. 



N- = C^e^^'^Sl. (38) 



using Eq. ()23|) we can calculate the corresponding transverse extrinsic curvature, which in 
the present case takes the form 

Kb = diag. |-C p gy ^ C— , C— sm 0| . (41) 
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Then, from Eqs. ()41|) and ()34p . we find the surface energy- momentum tensor ()26|) on the null 
hypersurface u = Uq, composed of the surface energy density a of the null shell and the 
pressures P in the 9- and 99-directions (J2Z|), as 



a 



Jw-h e[r)-m[r) 

V2(^ = -^aC -. — 2 , 

svrr 47rr"^ 

in 



in 



fl dV 



dr 



dm{r) 
dr 2a^(j)^ ' 40^02 



m[r)r r 



where m{r) is the mass function of the Vaidya region on the boundary u = Uq, 



(42) 



m{r) = m{V)\u= 



(43) 



Here we take the sign factor 772 to be 1 if the radiation is to the past of the wormhole region, 
and to be —1 if the radiation is to the future of the wormhole region (Fig|3J. In calculating 
the above equation, we have not used the particular expressions for the functions ip and /. 
Thus, it is valid generally in the case that the boundary surface is u =constant. 



Vaidya 




Vaidya 



Vaidya 




Vaidya 



FIG. 3: Matching domains: r/2 = 1 (left) and r/2 = — 1 (right). 



Henceforth we consider only the dust shell case P = 0, then we require 

dm mr r^ 

dr 2a?(j? Aa?(j? 



Integrating Eq. (j44j) . we obtain 



(44) 



m = ^(e-'' + 2/0(/) - 2<p{lfe'') + C0(/)e'', 
where C is an integration constant and related to a by 



(45) 



a = r]2CC 



Anr'^ 



(46) 
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Then if there is no hght-hke shell, a = and the mass function is continuous across the 
boundary surface, e = m. Extending the relation (jHBj) to the Vaidya region, introducing ziV) 
with z = I on the boundary surface, we obtain the mass function of the Vaidya solutions 
beyond the boundary surface 

m{z) = ^{e-'' + 2z(j){z) - 2<p{zfe'^) + C0(z)e^' (47) 

where the relation between z and V is 

v.,.) ^ -c /' ^"^--'t;-;' +^^^'^') .,. (48) 

j a(p{y) - C 

Transforming the energy-momentum tensor of the impulsive radiation, 

du Anr 



T- = ~t"'6{u - uo) = -t'-'6{V -Vo)— = ~V2CC^6{V - Vo). (49) 



C. Combined matching between Schwarzschild and static wormhole via Vaidya 

sSchwarzschild 
Vaidya \/ Vaidya Vaidya \( Vaidya 





FIG. 4: Matching domains: rji = r]2 = I (left) and ryi = r/2 = — 1 (right). 

In this subsection, we consider a collision of two oppositely moving impulses with Vaidya 
regions in opposite quadrants and Schwarzschild and static-wormhole regions in the other 
two quadrants. Connecting the impulses by (j3ip and (j49j) . we find that rji = ri2 = 1 ii the 
future region is wormhole, and rji = t]2 = —1 if the future region is Schwarzschild (Fig^. 
That is, rji = ri2 in both cases. Then the constant C is determined as 

vA 
Here we implicitly use the fact that the jump of a j ump vanishes, i.e. the jump across one 
impulse does not jump across the other impulse J2J]. 
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V=Vo 



FIG. 5: (Color online). Penrose diagram of the wornihole construction model. The wavy blue bold 
lines represent impulsive radiation with negative energy density. The region S is Schwarzschild, V 
is Vaidya and W is static-wormhole. The boundary between S and V corresponds to 2; = and 
the boundary between V and W corresponds to z = I. 

IV. WORMHOLE CONSTRUCTION FROM SCHWARZSCHILD BLACK HOLE 

In this section we present analytic solutions which describe the construction of a static 
wornihole from a Schwarzschild black hole. The whole picture is represented by Fig. El and 
the strategy is as follows. Firstly, impulsive phantom radiation is beamed in, causing the 
trapping horizons to jump inward, much as a shell of normal matter makes a black-hole 
trapping horizon jump outward. By controlling the energy and timing of the impulses, 
the trapping horizons can be made to instantaneously coincide. They can then form the 
throat of a static wornihole if constant-profile streams of phantom radiation are beamed in 
subsequently, with the energy density appropriate to a wormhole of that area. 

A. Vaidya region V 



First, we set up an initial Schwarzschild region S (jl3|) with mass M. Now we beam in 
impulsive phantom radiation symmetrically from either side, with the mass-energy of the 
shell being /i = 47rr^cr, then turn on constant streams of phantom radiation immediately 
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after the impulses. Then the region V should be Vaidya ()16|) with some mass function m, 
which on the boundary V = Vq between S and V is 

m{Vo) = M -r]Cfi = M + n, (51) 

from the matching formula between Schwarzschild and Vaidya ((221) • Here we must take 
(^ = — 1 and rj = 1, since the impulse is ingoing into the black hole, and the radiation is the 
future of the Schwarzschild region. 

In order for the final region W to be a static-wormhole region, the mass function m of 
the Vaidya region V must take the following form; 

m{z) = ^(e-^' + 2z(j){z) - 2<p{zfe'') + C<P{z)e'\ (52) 

and the relation between the coordinates V and z is 

From Eqs. (jM]) and (j3^ . a is 

a = 2(M + /i). (54) 

Connecting the impulsive radiation (j^ and (|mj) at / = 0, the constant C is decided as 

C ^ -^. (55) 

B. Wormhole region W 

We consider the spacetime in the region W. Using the matching formula ^^ between 
Vaidya p6|) with the mass function ()52|) and a wormhole, we find that the region W is a 
static wormhole with the mass-energy 

e{l) = {M + ij){e-'' + 2l(f){l) - 2<p{lfe''). (56) 

The relation between the throat radius tq of the wormhole in the final region W and the 
Schwarzschild mass M in the region S is 

ro = a = 2M + 2/x. (57) 
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Now we consider the energy and timing of the impulse. The tortoise coordinate r* inside 
a Schwarzschild black hole can be defined as 

r* = -r-2M\n(l- ^] for r < 2M, (58) 

so that dr*/dr > 0. In addition, the symmetry of the impulses means that the intersection 
point O is given by t = 0, r = ro or r* = Tq. Then the Eddington-Finkelstein relation (|Tl|l 
at the point O where the impulses collide gives 

Vo = -r* = 2(M + /i) + 2M In (^] . (59) 

Thus the energy and timing of the impulses are related. From this relation, first, the energy 
of the impulses must be always negative, /i < 0. The throat radius of the final wormhole 
(I57|) must be less than the horizon radius of the initial Schwarzschild black hole. Second, 
the later the negative-energy impulses occur, the larger the absolute value of the energy of 
the impulses must be. These features are consistent with the results of the 2D model 15|. 
In order for the final state not to have a naked singularity but to be a wormhole, the 
inequality —M < /i < is required. The throat radius of the wormhole in the final region W 
must be smaller than the horizon radius of the initial black hole, tq < 2M, since /i must be 
negative. In summary, one can prescribe the initial black-hole mass M > and the impulse 
energy /i G (— Af, 0) as free parameters, then the timing Vq of the impulses and the throat 
radius tq of the final wormhole are determined. 

V. WORMHOLE ENLARGEMENT BY IMPULSIVE RADIATION 

In this section we present an analytic solution which represents the enlargement of a 
static wormhole. The whole picture is represented by Fig. El and the strategy is as follows. 
Basically we want to open then close an expanding region of past trapped surfaces, by 
moving apart then rejoining the two trapping horizons comprising the wormhole throat in 
Wl. The general recipe is to first strengthen then weaken the negative energy density J8[. 
This can be done with a two-shot combination of primary impulses with negative energy, 
followed by secondary impulses with positive energy. To make the situation analytically 
tractable, the constant-profile phantom radiation is turned off between the impulses, leaving 
the region S as Schwarzschild and the regions VI, V2 as Vaidya. By controlling the energy 
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V=V2 



V=Vi 



FIG. 6: (Color online). Penrose diagram of the wormhole enlargement model. The wavy blue and 
red bold lines represent impulsive radiation with negative and positive energy density, respectively. 
The regions Wl and W2 are static-wormhole, VI and V2 are Vaidya, and S is Schwarzschild. The 
boundary between Wl and VI corresponds to z = I, the boundary between VI and S corresponds 
to z = 0, the boundary between S and V2 corresponds to w = 0, and the boundary between V2 
and W2 corresponds to w = I. 

and timing of the impulses and the energy density of the final constant-profile radiation, the 
final region W2 is also a static wormhole, but larger. 

A. Vaidya region VI 



We set up the initial region Wl as a static wormhole (jH)) with throat radius ri. Then the 
gravitational energy ei is 

2/0(/)-20(/)V'). 



ei(0 



— (e 
2 ^ 



(60) 



We beam in primary impulses symmetrically from both universes, then turn off the constant 
ghost radiation immediately after the impulses. Then the region VI should be Vaidya. 
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Timing the impulses at V^ = Vi, the matching formula ()46p between static- wormhole and 
Vaidya regions yields the mass-energy mi in the region VI as 

mi{z) = ^(e-"' + 2#(2) - 20(2)2e^') + Ci0(2)e^' (61) 

where the relation between the coordinates V and z is given by 

.. ,,;,-„■( .^, ,,(,)) ^^^^ (62) 

Jo nnv) - Cx 

Connecting the impulsive radiation from the boundary between W and VI to that between 
VI and S, we can decide the constant C\ = /ii/-\/A, where /ii is the mass-energy of the 
primary impulses. Then the mass function of the first Vaidya region VI at the boundary 
with Schwarzschild S becomes 

^i(^ = 0) = | (63) 

since the boundary surface V = V\ coincides with z = 0. 



B. Schwarzschild region S 

The region S is vacuum and therefore Schwarzschild. From the matching formula between 
Vaidya and Schwarzschild ()29p . the mass M of the Schwarzschild region S becomes 

M=^ + r7iCi/ii = y-Ci/ii, (64) 

where 771 = —1, since the Vaidya region VI is to the past of the Schwarzschild region S. Since 
we construct solvable symmetric models, the impulses are also both ingoing or outgoing. This 
means the region S must be inside a black-hole or white-hole region and r\ < 2M. So when 
the impulse has negative energy, /ii < 0, the sign of Ci must be positive, 

M='^-^^,. (65) 

This means the radiation is outgoing, 

dr , . 

and S is a white-hole region. This is what we need to enlarge the wormhole, as a white-hole 
region is expanding. Conversely, one could reduce the wormhole size by taking ^i > 0, 
creating a contracting black-hole region. 
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Now the symmetry of the impulses means that the intersection point Oi is given by 
t = 0. Then the Eddington-Finkelstein relation (|T^ at the point Oi where the impulses 
collide gives 

V, = -rl = n + 2Mln (l - ^) = 2(M + /xi) + 2Mln (^^^ . (67) 

Eqs. ()64p and ()67|) mean that the timing Vq of the impulses and the Schwarzschild mass 
M are determined by the throat radius ri of the initial wormhole and the energy /ii of the 
impulses. 

C. Vaidya region V2 

We next beam in secondary impulses symmetrically from both universes, and turn on 
constant-profile phantom radiation immediately after the impulses. Then the region V2 
must be Vaidya. Timing the impulses at V^ = V2, the matching formula between Vaidya and 
Schwarzschild fj29|) yields the mass function of the second Vaidya region V2 as 

m2{V2) = M - 7]2C2fi2 (68) 

on the boundary, where fi2 is the mass of the second shell. The region V2 is an outgoing 
Vaidya region which is to the future of the Schwarzschild region, so that ri2 = I and ^2 = 1- 
Then 

^2(^/2) = Af-/i2 = y-/ii-/i2. (69) 

Here in order for the final region W2 to be a static wormhole, the mass function 7712 must 
take the following form, 

m2{w) = ^(n - 2/ii - 2/i2)(e-"'' + 2i/;0(«;) - 20(^)^6"'') + C20(w)e"'', (70) 

where the coordinate w is related with V by 

r 2rie->Y^2ym) ,^^y^^ (71) 

Jo r2(f){y) - C2 

from the matching formula ()45p . Here r2 is the throat radius of the wormhole in the final 
region W2. Connecting the impulsive radiation from the boundary between the regions S 
and V2 to that between V2 and W2, we can decide the constant C2 = /i2/ vA- 
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It can be shown that there are trapping horizons in the regions V2, as depicted in Fig 
El the negative-energy and positive-energy impulses respectively make the horizons jump 
to the future and the past. It is difficult to study the horizons analytically in the Vaidya 
coordinates, but it can be shown that they are null using dual-null coordinates x^, as follows. 
If we have V pointing along x~^, then there is only a T++ component in the energy tensor. 

The T and T+_ components of the Einstein equations [lO| then show respectively that, 

where 9_r = 0, then d^d^r = and d^d^r < 0, which means that the horizon d^r = is 
null. Similar behavior occurs in the 2D model [l5[ , though the horizons were omitted in the 
corresponding diagram. 

D. Wormhole region W2 

Finally, we consider the spacetime in the region W2. Since constant-profile phantom 
radiation is beamed in for \^ > V2 in order for the region V2 to be Vaidya (fTT)|) with the 
mass function ()7()j) . we can match it to a static- wormhole region from the matching formula 
PH|l . We find that the mass function €2 of the wormhole in the region W2 is 

62(1) = in - 2^1 - 2fi2){e~'" + 2/0(/) - 20(/)2e''), (72) 

and the throat radius r2 is 

rs = ri - 2/ii - 2/i2 = 2M - 2/22- (73) 

Again, the symmetry of the impulses means that the intersection point O2 is given by 
t = 0. Then the Eddington-Finkelstein relation (fTHl at the point O2 where the impulses 
collide gives 

V2 = -r* =r2 + 2Mln (l - ^) = 2(M + /i2) + 2Mln (g) . (74) 

From this relation, the energy of the impulses must be positive, ^2 > 0. In addition, the 
inequality 

r2 > ri (75) 

holds, since the region S is part of a white-hole region. That is, the wormhole is enlarged. 
We find that the absolute value |yUi| of the energy density of the primary impulses should 
be stronger than that of the secondary impulses, 

lAiil > 1/^2 1 (76) 
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from Eqs. d7H|) and d7H|l . 

We find the relation between the energy and timing of impulses as 



^2 - ^1 = rl - r{ = 2{fii + /is) + 2M In (^j , (77) 



from Eqs. ()fi7|l and ()74j) . Eq. ()77|1 means that the longer the interval between the first and 
second impulses is, the smaller the value of energy density of the second impulse must be. 
In summary, once the throat radius of the initial wormhole and the energy of the impulses 
(ri,yUi,/i2) are prescribed, the timings of the impulses, the intermediate Schwarzschild mass 
and the throat radius of the final wormhole (Vi, M, J^rs) are determined. These features 
are also consistent with the results of the 2D model 






VI. JUMP IN ENERGY DUE TO IMPULSIVE RADIATION 

A general spherically symmetric metric can be written in dual-null form as 

ds"^ = r^dn"^ - h dx+dx- (78) 

where r > and h > are functions of the future-pointing null coordinates (x"*", x~). Writing 
d± = d/dx^, the propagation equations for the energy E (0) are obtained from the Einstein 
equations as |uj 

d±E = 87rh-^r\T+_d±r - T±±d^r). (79) 

We have considered impulsive radiation defined by 

Ta. = ^^^(-" - -o) (80) 

where the constant Xq gives the location of the impulse and the constant /i± is its energy. 
More invariantly, the vector p = ~g^^{fJ'±dx^) is the energy-momentum of the impulse. 
Then the jump 

pXQ+a 

[E]± = lim / d±Edx^ (81) 



°-0^xo-a 



in energy across the impulse is given by the jump formula 

[E]± = c^fi±, c± = -2ih-^d^r)\,±=,,. (82) 

The vector c = c^9+ + c^9__ is actually c = g^^{dr) and so 

[E]± = -p-dr (83) 
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is a manifestly invariant form of the jump formula. Note that while the energy-momentum 
vector p (or ^±dx^) is invariant, the energy /i± depends on the choice of null coordinate 
x^, reflecting the fact that a particle moving at light-speed has no rest frame and no pre- 
ferred energy. However, in a curved but stationary space-time, the stationary Killing vector 
provides a preferred frame and a preferred energy /i±. 

We need only employ the jump formula in the following cases, (i) Inside a Schwarzschild 
black hole, we can take future-pointing x^ = — r* ± t where dr* jdr = (2M/r — 1)"^. Then 
r* = _(a;+ + x")/2, h = 2M/r-l and9±r = {dr/dr*)d±r* = (l-2M/r)/2 gives c=^ = 1 and 
[E]± = fi±. (ii) Inside a Schwarzschild white hole, we can take future-pointing x^ = r* ±t 
and similarly obtain c^ = — 1 and [E]± = —fi±. (iii) On the throat of a static wormhole, 
where d+r = d^r = and h is finite iJj, one finds c^ = and [E]± = 0. This is summarized 

as 

I/i inside a Schwarzschild black hole 

— yU inside a Schwarzschild white hole (84) 

on the throat of a static wormhole 

where the indices on fj,± are now omitted. 




FIG. 7: an infinitesimal box across which two radiative impulses (arrows) intersect. 

Now assuming an infinitesimal diamond-shaped box around the point where the impulses 
collide as in FigQ we can evaluate the jump in energy across the impulses for each case 
in Sees. IV and V. For wormhole construction (FiglSj), the energy E will jump by fi from 
the region S to V and by from the region V to W, evaluated in the limit at the point, 
recovering ro/2 = M + ^ ()57|) . where M is the black-hole mass in the region S. Similarly 
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for wormhole enlargement (FiglHI), the energy E will jump by from Wl to VI and by —^i 
from VI to S at the point Oi, and by — /i2 from the region S to V2 and by from the region 
V2 to W2 at O2, recovering M = ri/2 - fn §B> and r2/2 = M - ^2 dZSI), where M is the 
black- hole mass in the region S. 

Thus even without performing the detailed matching, the basic properties of the solutions 
could be predicted simply by the jump formula for E and continuity of the area A = Anr"^. 
For wormhole construction, continuity of r at O implies Tq < 2M, and the jump formula 
gives /i < 0; the impulses must have negative energy. Similarly, for wormhole enlargement, 
continuity of r at 01 and 02 implies ri < 2M and r2 < 2M, and the jump formula gives 
/ii < and yU2 > 0; the primary and secondary impulses must have negative and positive 
energy respectively. 

VII. SUMMARY 

In this paper, we have studied wormhole dynamics in Einstein gravity under (phantom 
and normal, impulsive and regular) pure radiation, constructing analytic solutions where 
a traversable wormhole is created from a black hole, or the throat area of a traversable 
wormhole is enlarged or reduced, the size being controlled by the energy and timing of the 
impulses. The solutions are composed of Schwarzschild, static-wormhole [lJ| and Vaidya 
regions matched across null boundaries according to the Barrabes- Israel formalism. For this 
purpose we have derived the matching formulas which apply when the direction of radiation 
in the Vaidya region is either parallel or transverse to the boundary. These formulas are 
useful for other problems. 

The results provide concrete examples of how to create and enlarge traversable wormholes, 
given the existence of Schwarzschild black holes and phantom energy. We have worked within 
standard General Relativity, inventing no new theoretical physics other than an idealized 
model of phantom energy, on which general arguments do not depend j^ . Thus if space-time 
foam and phantom energy do exist and can be controlled, then traversable wormholes can 
be constructed and enlarged. 
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